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Abstract 

Starting from basic physical principles, we present a novel derivation 
linking the global wind power to measurable atmospheric parameters. The 
resulting expression distinguishes three components of the atmospheric po¬ 
wer (the kinetic power associated with horizontal and vertical motion and the 
gravitational power of precipitation) and highlights problems with previous 
approaches. Focusing on Laliberte et al. (2015), we show how inappropriate 
treatment of material derivatives in the presence of phase transitions leads 
to significant errors in wind power analyses. We discuss the physical con¬ 
straints on global wind power and the opportunities provided by considering 
the dynamic effects of water vapor condensation. 

1 Introduction 

Atmospheric scientists have long sought fundamental principles which can ex¬ 
plain global circulation power. This search has gained renewed significance in 
view of the need to reconcile an apparent mismatch between model predictions 
and observed trends (e.g., Kociuba and Power, 2015). While models suggest cir¬ 
culation should slow as global temperatures increase, independent observations 
indicate that in fact global wind power is increasing (e.g., de Boisseson et ah, 
2014). Without a coherent theoretical framework such results cannot be readily 
reconciled. Meanwhile no expression for global atmospheric power in a moist 
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atmosphere has been derived from fundamental physical principles. Here we 
derive such an expression. This permits us to examine the determinants of the 
atmospheric power budget and to identify and address some widespread misun¬ 
derstandings. In particular we argue that equilibrium thermodynamics offer little 
insight into the constraints on global circulation power. 


2 What is atmospheric power? 


2.1 Total power 

Work in the atmosphere is done by expanding air. Atmospheric water in its solid 
and liquid form is incompressible and neither performs work nor occupies any 
appreciable volume (Pelkowski and Frisius, 2011). Work per unit time (power) of 
an air parcel containing N moles and occupying volume V (m^) is 


dV ~dV ^dN 


dp n \ 


dN 


pV— = N -V-f- + R— + RT—. 


dt 


dt J 


dt 


( 1 ) 


Here we have used the ideal gas law, pV = RT and pdV = RdT — Vdp, where 
T is temperature, V = N~^ is the atmospheric volume occupied by one mole of 
air, N is air molar density (mol m“^), V = NV, p is air pressure and R = 8.3 
J mol“^ is the universal gas constant. 

In Eq. ([T]) the expression in braces represents part of the parcel’s power per 
mole that is independent of the rate of phase transitions dN/dt. To express it per 
unit volume we use the relationship N = VN. Then total power Wtot (W) in 
an atmosphere containing n air parcels with total fixed volume V = ^ = 

fy, dV (subscript i refers to the i-th air parcel) is 


Wtot = 


^ dVi f / dp ^^dT -A 

> pi^ = / —^ + NR— + RTN 

^ dt /v \ dt dt / 


dV. 


( 2 ) 


Here N is the molar rate of phase transitions per unit volume (mol m“^ s“A- Its 
integral over volume V is equal to the total rate of phase transitions in all the n air 
parcels: NdV = dNi/dt. 

We will now use the definition of material derivative for X = {p,T, N}: 


the continuity equation: 


dX dX _ 

(3) 

dN 

V.(iVv) = iV-^. 

(4) 
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and the divergence theorem: 


/ NY-VTdV= I NT{Y-n)dS- [ TV-(iVv)dV=- / TV-(iVv)dV. (5) 
Jv Js Jv Jv 

Here v is air velocity, n is unit vector perpendicular to the unit surface of area diS; 
the integral is taken over the surface enclosing the atmospheric volume. Since the 
air does not leave the atmosphere, we have v ■ n = 0, and the integral over S in 
([5]) is zero. With help of Eqs. Q-® and the ideal gas law dp/dt = Rd{NT)/dt 
we obtain from Eq. Q 

Wtot = - / (v ■ Vp)dV = - / (u ■ Vp + w • Vp)dV, (6) 

Jv Jv 

where u and w are the horizontal and vertical air velocities, v = u + w. 

2.2 Kinetic power 

In Eq. ® the expression — (v ■ Vp) represents work performed per unit time per 
unit air volume by the pressure gradient. The horizontal pressure gradient gen¬ 
erates the kinetic energy of the horizontal wind. The vertical pressure gradient 
generates the kinetic energy of the vertical wind plus it changes the potential en¬ 
ergy of air in the gravitational field. 

In hydrostatic equilibrium of air we have 


VzP = pg, (7) 

where p = NM is air density (kg m“^), M is air molar mass (kg mol“^). In 
the real atmosphere due to the presence of non-gaseous water the air distribution 
deviates from Eq. © such that we have V zP = (p+Pz)g and in Eq. ® — w ■ Vp = 
—pw ■ g — p;w ■ g, where pi is mass density of the non-gaseous water in the air. 

Term — pw • g represents the vertical flux of air: it is positive (negative) for the 
ascending (descending) air. Recalling that g = —gVz and using the divergence 
theorem and the stationary continuity equation, see @, we can write 

Wp = — pw-gdV = / n-{pwgz)dS— / gzV■{pw)dV = — gzpdV. (8) 
Jv Js Jv Jv 

Eor a dry atmosphere where p = 0, this last integral in Eq. ([8]) is zero and Wp = 0: 
indeed, in this case at any height 2 ; there is as much air going upwards as there is 
going downwards. In a moist atmosphere, evaporation p > 0 makes a negative 
contribution to Wp, while condensation p < 0 makes a positive contribution. 
This is because Wp reflects the work of water vapor as it travels from the level 
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where evaporation oeeurs (where water vapor arises) to the level where eonden- 
sation oeeurs (where water vapor disappears). When eondensation oeeurs above 
where evaporation occurs, the water vapor expands as it moves upwards towards 
condensation, and the work is positive. If condensation occurs below where evap¬ 
oration occurs, the water vapor must compress to reach the evaporation site; thus, 
the work is negative. When evaporation occurs at the Earth’s surface z = 0, Wp 
dS]) is equal to PgHpSp, where Hp is the mean height of condensation, P is 
global mean precipitation at the surface (kg m~^ s“^) and Sp is Earth’s surface 
area (Makarieva et ah, 2013a). It is natural to call Wp the ’’gravitational power of 
precipitation”. 


Thus for the stationary total power we have 



Wtot 

= Wk + Wp, 


( 9 ) 

Wk 

= - (u ■ Vp-f pzw • g)dV 

Ri — u - VpdV, 

(10) 


Jv 

Jv 


Wp 

= - [ gzpdV ^ PgHpSp, 

p = - [ pdV/Sp. 

(11) 


Jv 

Jz>0 


Per unit area, global Wp and Wp were estimated at 2.5 and 0.8 W m 

respec- 


tively (Huang and McElroy, 2015; Makarieva et ah, 2013a). 

The term —p/w • g in Eq. (fTOl) is not related to the gravitational power of 
precipitation. It describes kinetic energy generation on the vertical scale of the 
order of the atmospheric scale height H = —p/{dp/dz) = RT / {Mg) "zz 10 km. 
This energy is generated because the vertical air distribution deviates from the 
hydrostatic equilibrium ©. Hydrometeors act as resistance not allowing the non¬ 
equilibrium pressure difference Ap ~ PigH to be converted to the kinetic energy 
of the vertical wind. In the atmosphere on average pi/p ^ 10“^ (Makarieva et ah, 
2013a). Without hydrometeors, a pressure difference Ap ~ IW^pgPl ~ 1 hPa 
would produce a vertical velocity of about te ~ 1 m s“^ (pu'^/2 = Ap). This 
is two orders of magnitude larger than the characteristic vertical velocities w ~ 
10“^ m s“^ of large-scale air motions. Hydrometeors thus have an effect similar 
to turbulent friction at the surface which does not allow horizontal velocities to 
develop. Eor example, the observed meridional surface pressure differences of the 
order of Ap ~ 10 hPa in the tropics, if friction were absent, would have produced 
horizontal air velocities of about 40 m s“^ instead of the observed 7 m s“^ The 
term —p/w ■ g is less than 1% of Wtot and can be neglected: its volume integral 
taken per unit surface area is less than pigHw ~ IW^pw ~ 10“^ W m“^, where 
p = pgH = 10^ Pa is air pressure at the surface. 
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3 Revisiting the current understanding of the atmo¬ 
spheric power budget 

3.1 The physical meaning of Wtou and Wp 

To our knowledge, Gorshkov (1982) was the first to estimate Wp for land assum¬ 
ing Hp = 2 km. In the meteorologieal literature, Pauluis et al. (2000) defined 
preeipitation-related frietional dissipation as Wp + Jy piw ■ gdV and estimated its 
value for the tropies. This estimate was later revised by Makarieva et al. (2013a). 
Pauluis et al. (2000) stated that in the stationary ease total power of atmospheric 
motions must be equal to the sum of the dissipation rate at the microscale around 
the hydrometeors and the rate of dissipation of motions on convective scale. The 
latter was approximated it as the power of the buoyancy force that is proportional 
to vertical velocity (Pauluis et al., 2000, their Eq. 8). Derivation of Eqs. (l9l)-(fTTI) 
have not been previously published. 

These equations provide insight into the controls on circulation power. Eirst, 
we can see that Wp does not depend on the interaction between the air and the 
falling drops. This term would be present in the atmospheric power budget even 
if drops were experiencing free fall and did not interact with the air at all (such 
that no frictional dissipation on drops occurred in the atmosphere). Second, the 
small term —piw ■ g describes generation of kinetic energy at the convective scale 
H (while this energy dissipates on hydrometeors) and thus should be included 
into the total rate of kinetic energy generation Wk- Einally, in the absence of a 
derivation relating Eq. ([T]) and Eq. ([91), it is not clear a priori whether precipitation 
makes an independent contribution to the total power budget. Indeed it could 
appear as a form of turbulent dissipation of the kinetic energy generated at larger 
scales (i.e. water is lifted by turbulent diffusion), and thus be comprised in the 
expression for the total kinetic power. 

We disagree with the recent commentary of Pauluis (2015) on the work of 
Ealiberte et al. (2015). To provide some context, an ideal atmospheric Carnot 
cycle consuming heat flux F = 100 W m“^ {F is limited from above by solar 
power flux reaching the planetary surface) at surface temperature Tin = 300 K 
and releasing heat at Tout = Tin — AT with AT = 30 K being roughly the 
magnitude of the Earth’s greenhouse effect, would generate kinetic energy at a 
rate of Wc = T(AT/Tj„) = 10 W m“^. Ealiberte et al. (2015) estimated total 
atmospheric power Wtot per unit area to be around 4 W m~^. Comparing this 
result with Wc, Pauluis (2015) notes that ’’estimates for the rate of kinetic energy 
production by atmospheric motions are about half this figure”. Confusion has 
apparently arisen between total atmospheric power Wtot and kinetic power Wk- 
Indeed, Pauluis (2015) continues that ’’the difference is very likely due to Earth’s 
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hydrological cycle, which reduces the production of kinetic energy in two ways”, 
one of whieh is the gravitational power of preeipitation Wp. However, as is elear 
from Eq. @, Wp is equal to the difference between Wtot and Wk and is unrelated 
to the differenee Wc — Wtot- 

Since the relationship between Wc, Wp and Wk was discussed by Pauluis 
et al. (2000) (see their Eq. (7)), the confusion may stem from a misinterpreta¬ 
tion of the magnitude estimated by Laliberte et al. (2015) as atmospherie power, 
Jj^{l/p){dp/dt)dM = Jy{dp/dt)dV, for Wk - while it was in faet Wtot- No¬ 
tably, neither Pauluis et al. (2000) nor Ealiberte et al. (2015) identified the faet 
that Wk depends on horizontal and not vertieal veloeities. This is essential for 
eomparing theory and observations. Horizontal pressure gradients and wind ve¬ 
loeities are observed, while the vertieal veloeities are inferred with signifieant 
uneertainty. 

There are other studies where kinetie energy generation is estimated from hor¬ 
izontal veloeities as Wk flO]) (see, e.g., Boville and Bretherton, 2003; Huang 
and McElroy, 2015). At the same time, Wk is sometimes confused for the to¬ 
tal atmospheric power: i.e. in the total power budget the gravitational power of 
precipitation, Wp, is overlooked (e.g., Huang and MeElroy, 2015, their Eig. 10). 


3.2 The analysis of Laliberte et al (2015) 


Equations (fTh- tfTT]) also allow us to identify errors in analyses of Wtot- To estimate 
the atmospherie power budget, Ealiberte et al. (2015) adopted the thermodynamie 
identity 


T— = 


dh 

dt 


1 dp 
p dt 


p- 


dqp 


( 12 ) 


dt dt p dt dt 
where s is entropy, h is enthalpy, p is ehemieal potential (all per unit mass of wet 
air), 1/p is specifie air volume and qp is water mixing ratiqli When integrating 
Eq. (fT^ over atmospherie mass, Laliberte et al. (2015) noted that the enthalpy 
term vanishes, J^{dh/dt)dM = 0, beeause the atmosphere is in a steady state. 
However, Eqs. O]) and (HI) indieate that for any sealar quantity, in our ease enthalpy 
h, in the stationary ease, dh/dt = 0, we have 



Here p is the mass rate of phase transitions per unit volume (kg m“^ s“^), V ■ 
(pv) = p. This integral is not zero in an atmosphere where phase transitions take 
plaee. 

'The unconventional sign at the chemical potential term follows from p being defined in 
Eq. (flD i relative to dry air: hence, when the dry air content diminishes (while dqr/dt > 0) 
this term is negative. For details see p. 8 in the Supplementary Material of Laliberte et al. (2015). 
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This omission creates a significant error. To obtain an approximate estimate 
we can assume that all evaporation occurs at the surface, while all condensation 
occurs at a global mean height of Hp = 2.5 km, where about half of surface 
water vapor has condensed in the ascending air (see Makarieva et ai, 2013a). 
Then we have [Jj^{dh/dt)dM]/SE ~ —P[h{0) — h{Hp)], see Eq. (fTTI) . where 
h{0) and h{Hp) are mean enthalpies at the surface {z = 0) and at mean con¬ 
densation height {z = Hp). With h ^ CpT + Lq, where q is water vapor mass 
mixing ratio, Cp = 10^ J kg“^ is heat capacity of air at constant pressure, 
L = 2.5 X 10® J kg“^ is latent heat of vaporization, T(0) — T{Hp) 15 K, 
g(0) — q{Hp) = 0.5g(0) ~ 0.5 x 10“^, P = 10^ kg m~^ yr“\ we find that 
lf^{dh/dt)dM]/SE rsj — 1 W m This figure is about one quarter of the total 
atmospheric power Wtot/Sp ~ 4 W m“^ estimated by Laliberte et al. (2015) for 
the MERRA re-analysis and the CESM model. 

As Eq. (fT^ is an identity, whether the incorrect assumption has increased or 
decreased the value of atmospheric power depends on the particular method of 
calculating it. Ealiberte et al. (2015) first calculated the mass integral of Tds/dt 
from Eq. (fT^ . then calculated fidqp/dt from atmospheric parameters and then 
used the obtained values to estimate the total power — Jj^{l/p){dp/dt)dAi. In 
such a procedure, putting Jj^{dh/dt)dM. = 0 resulted in an overestimate of the 
atmospheric power by about 1 W m“^. Since the omitted term depends on pre¬ 
cipitation rate, its omission is crucial not only for a correct estimate of the mean 
value of Wtot, but also for the determination of any trends related to precipitation, 
which is a focus of Ealiberte et al. (2015). Thus the quantitative conclusions of 
Ealiberte et al. (2015) appear invalid. 


4 Dynamic versus thermodynamic approach 

Under an assumption of local thermodynamic equilibrium Eq. (fT^ defines ds/dt 
via measurable atmospheric variables. As such, it does not carry any additional 
information about the atmosphere besides that contained in the set of local val¬ 
ues of pressure, temperature, relative humidity and air velocity. Total atmospheric 
power can be estimated from Eqs. (l^-(fTTI) without involving entropy. This fact il¬ 
lustrates a limitation of the thermodynamic approach in that it cannot explain why 
a particular thermodynamic cycle exists: i.e. why does the atmosphere generate 
power. 

Indeed, the kinetic power Wk can be viewed as a measure of the dynamic 
disequilibrium of the Earth’s atmosphere. In equilibrium, for example, under 
conditions of hydrostatic and geostrophic or cyclostrophic balance, no power is 
generated and Wk = 0. Considerations of entropy allow one to constrain the 
maximum work that can be extracted from a given thermodynamic cycle. How- 
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ever, equilibrium thermodynamics cannot predict whether the considered system 
will be in equilibrium. It thus provides no insight as to if or why Wk on Earth 
differs from zero. The obvious upper limit - the global efficiency of solar energy 
conversion into useful work, which amounts to about 90% (Wu and Liu, 2010; 
see Pelkowski, 2012 for a rigorous theoretical discussion) - appears irrelevant for 
constraining Wk'- Wk is far below this thermodynamic limit, being about 1% of 
incoming solar radiation. 

In contrast, the dynamic approach identifies situations when production of 
kinetic energy cannot be zero in principle. Saturated water vapor in the grav¬ 
itational field above the liquid ocean is a unique dynamic system, because its 
pressure is controlled by temperature alone rather than by temperature and mo¬ 
lar density as for the non-condensable ideal gases. Specifically, while dry air 
can rise in hydrostatic equilibrium, the saturated water vapor cannot. In an at¬ 
mosphere composed of pure saturated water vapor the water vapor rises (return¬ 
ing to the Earth as precipitation) and kinetic energy is produced at a local rate 
—w{dp.u/dz — p-uq) > 0, where is mass density of water vapor (Makarieva 
et ah, 2015b). If non-condensable gases are added to such an atmosphere, it be¬ 
comes possible to arrange a hydrostatic equilibrium in the vertical plane, such 
that all the kinetic power that derives from condensation is now generated in the 
horizontal plane: —u ■ Vp = —w{dpy/dz — p^/H) = —wRTNd'y/dz. Since 
— J^^QwN{d'y/dz)dz ~ P/My, where My is molar mass of water vapor, the 
global kinetic power that can be derived from condensation is about 4 W m“^ (as¬ 
suming mean temperature of condensation T 270 K and P = 10^ kg m-^ yr-^) 
(Makarieva et al, 2013a). This theoretical estimate is about 40% higher than 
the most recent estimate of Wk/Se = 2.5 W m~^ obtained using the MERRA 
re-analysis (Huang and McElroy, 2015). 

The difference is likely related to an insufficient spatial and temporal resolu¬ 
tion. Even when in real eddies there is no evaporation in the descending air (i.e. 
for te < 0 we have d'^/dz = 0 and N = 0), since in the rising saturated air we 
always have d'-y/dz < 0, the spatial and temporal averaging produces a non-zero 
vertical gradient d'y/dz < 0 everywhere. This results in ’’spurious” evaporation 
N Ri wd'y/dz > 0 for z > 0. It diminishes the integral of — J^^^NdV and 
all the related quantities, including the kinetic power derived from condensatior@. 
Our theoretical estimate Wk = 4 W m“^ suggests that with increasing spa¬ 
tial and temporal resolution the atmospheric power estimated from observations 
should grow. In agreement with these ideas, Kim and Kim (2013) using the daily 
averaged MERRA data obtained Wk/Se = 2 W m“^, which is 20% less than the 
value obtained by Huang and McElroy (2015) using 3-hour observations (Huang 

^This spurious evaporation also diminishes the integral of dh/dt (fT3l l. a problem not discussed 
by Laliberte et al. (2015). 
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J, personal communication). 

In global circulation models a non-zero rate of kinetic energy generation is 
achieved by introducing an ad hoc intensity of turbulent diffusion which deter¬ 
mines the rate at which kinetic energy is dissipated (and, in the steady state, gen¬ 
erated) (see Makarieva et ah, 2015a). It is this parameterization that postulates 
a certain value of Wk and controls its behavior. Since this parameterization is 
unrelated to the hydrological cycle (i.e. one and the same turbulent diffusion co¬ 
efficient can be used in both dry and moist models), comparing Wk across models 
with varying intensity of the hydrological cycle does not shed light on the actual 
role of water vapor. Meanwhile the condensation-induced atmospheric dynamics 
outlined above (for further details see Makarieva et ah, 2013b; 2015b and ref¬ 
erences therein) suggests that in the absence of the hydrological cycle the atmo¬ 
spheric power on Earth would have been negligible compared to what it is today. 
We thus urge attention to the dynamic effects of condensation. 


Acknowledgements 

We thank Dr. Olivier Pauluis, Dr. Frederic Laliberte and five anonymous reviewers 
for their critical comments. This work is partially supported by Russian Scientific 
Foundation under Grant no. 14-22-00281 and the University of California Agri¬ 
cultural Experiment Station. 

References 

Boville BA, Bretherton CS. 2003. Heating and kinetic energy dissipation in the 
NCAR Community Atmosphere Model. Journal of Climate 16: 3877-3887. 
de Boisseson E, Balmaseda MA, Abdalla S, Kallen E, Janssen PAEM. 2014. How 
robust is the recent strengthening of the Tropical Pacific trade winds? Geophysical 
Research Letters 41: 4398-4405. 

Gorshkov V.G. 1982: Energetics of the biosphere. Eeningrad Politechnical Insti¬ 
tute, 80 pp. 

Huang J, McElroy MB. 2015. A 32-year perspective on the origin of wind energy 
in a warming climate. Renewable Energy 77: 482-492. 

Kim YH, Kim MK. 2013. Examination of the global lorenz energy cycle using 
MERRA and NCEP-reanalysis 2. Climate Dynamics 40: 1499-1513. 

Kociuba G, Power SB. 2015. Inability of CMIP5 models to simulate recent 
strengthening of the Walker circulation: implications for projections. Journal 
of Climate 2%\ 20-35. 

Ealiberte F, Zika J, Mudryk E, Kushner PJ, Kjellsson J, Doos K. 2015. Con¬ 
strained work output of the moist atmospheric heat engine in a warming climate. 
Science 347: 540-543. 


9 


Makarieva AM, Gorshkov VG, Nefiodov AV, Sheil D, Nobre AD, Bunyard P, 
Li BL. 2013a. The key physieal parameters governing frietional dissipation in a 
precipitating atmosphere. Journal of the Atmospheric Sciences 70: 2916-2929. 
Makarieva AM, Gorshkov VG, Sheil D, Nobre AD, Li BL. 2013b. Where do 
winds come from? A new theory on how water vapor condensation influences at¬ 
mospheric pressure and dynamics. Atmospheric Chemistry and Physics 13: 1039- 
1056. 

Makarieva AM, Gorshkov VG, Nefiodov AV, Sheil D, Nobre AD, Li BL. 2015a. 
Comment on ’’The tropospheric land-sea warming contrast as the driver of tropical 
sea level pressure changes”. Journal of Climate 2^: 4293-4307, doi: 10.1175/JCLI- 
D-14-00592.L 

Makarieva AM, Gorshkov VG, Nefiodov AV. 2015b. Empirical evidence for 
the condensational theory of hurricanes. Physics Letters A, 379: 2396-2398, 
doi: 10.1016/j .physleta.2015.07.042. 

Pauluis O, Balaji V, Held IM. 2000. Frictional dissipation in a precipitating atmo¬ 
sphere. Journal of the Atmospheric Sciences 57: 989-998. 

Pauluis OM. 2015. Atmospheric science. The global engine that could. Science 
347: 475-476. doi: 10.1126/science.aaa368L 

Pelkowski J, Frisius T. 2011. The theoretician’s clouds - heavier or lighter than 
air? On densities in atmospheric thermodynamics. Journal of the Atmospheric 
Sciences 68: 2430-2437. 

Pelkowski J. 2012. Of entropy production by radiative processes in a conceptual 
climate model. Meteorologische Zeitschrift 21: 439-457. 

Wu W, Fiu Y. 2010. Radiation entropy flux and entropy production of the Earth 
system. Reviews of Geophysics 48: RG2003. 


10 


